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ABSTRACT 

We have analyzed two epochs of HST/WFPC2 observations of the young Galactic starburst cluster 
in NGC 3603 with the aim to study its internal dynamics and stellar population. Relative proper 
motions measured over 10.15 yrs of more than 800 stars enable us to distinguish cluster members 
from field stars. The best-fitting isochrone yields ^4v=4.6-4.7mag, a distance of 6.6-6.9 kpc, and 
an age of lMyr for NGC 3603 Young Cluster (NYC). We identify pre- main-sequence/main-sequence 
transition stars located in the short-lived radiative-convective gap, which in the NYC occurs in the 
mass range 3.5-3.8 M Q . We also identify a sparse population of stars with an age of 4Myr, which 
appear to be the lower mass counterparts to previously discovered blue supergiants located in the 
giant H II region NGC 3603. For the first time, we are able to measure the internal velocity dispersion 
of a starburst cluster from 234 stars with I < 18.5 mag to <7 pm irj = 141±27/ias yr -1 (4.5 ±0.8 km s _1 
at a distance of 6.75 kpc). As stars with masses between 1.7 and 9M all exhibit the same velocity 
dispersion, the cluster stars have not yet reached equipartition of kinetic energy (i.e., the cluster is 
not in virial equilibrium). The results highlight the power of combining high-precision astrometry and 
photometry, and emphasize the role of NYC as a benchmark object for testing stellar evolution models 
and dynamical models for young clusters and as a template for extragalactic starburst clusters. 
Subject headings: astrometry - open clusters and associations: individual (NGC 3603) - stars: evolu- 
tion - stars: formation - stars: pre-main sequence - stars: kinematics and dynamics 



1. INTRODUCTION 

Massive young stellar clusters are outstanding ob- 
jects containing copious numbers of stars over the en- 
tire stellar m ass range. With m as ses between 10 4 M^ 
and 10 7 M^ (IZhang fc Fall 19991 : Ide Griis et al. 20031: 
iMengel et al. 20081 ) . thev cover the upper end of the 
cluster mass function and may constitute progenitors of 
globular clusters (GCs; Zhang & Fall 1999; McCrady 
& Graham 2007). While extragalactic starburst clus- 
ters, such as th ose in the Antennae Galaxies , are often 
barely resolved ([Whitmore fc Schwcizer 1995), in Milky 
Way starburst clusters thousands of individual stars can 
be observed. 

In addition to three clusters in the Galactic Cen- 
ter region (Arches, Quintuplet, Young Nuclear Clus- 
ter), only a handful of Milky Way starburst clusters 
located in spiral arm s have so far been identified (e.g. 
iBrandner et al. 2008f ). Among the spiral arm clus- 
ters, the NGC 3603 Young Cluster (NY C), located in it s 
namesake giant Hn region NGC 3603 (jKennicutt 1984( 1. 
is the most co mpact and youngest cluster with an 
age of » lMvr (iBrandl et al. 19991 iStolte et al. 20041 : 
iSung fc Bcsscl i""2"b04l) and a centra l density po > 6 • 
10 4 M o pc" 3 (jHaravama et al. 20081 ). It hosts three 
Wolf- Rayet stars, at leas t 6 02/03, and 30 late O-type 
stars (M offat et al. 2004T) . and is extensively referenced 
as a template for extra galactic starburst environments 
(e.g. iLamers et al. 2006f) . 

Previously, dynamical studies of Galactic starburst 
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clusters were largely restricted to one-dimensional ve- 
locity dispersions derived from radial velocity measure- 
ments of_a_Je_w^_of_tlie_jnostJjan^ cluster members 
(e.g. IMengel fc Tacconi-Garman 20071 ). Using multi- 
epoch observations of GCs, King & Anderson (2001, 
2002) and Anderson & King (2003a) pioneered high pre- 
cision proper motion studies with the Hubble Space Tele- 
scope (HST), which enabled the distinction of GC mem- 
bers from field stars, and the study of GC dynamics and 
kinematics. 

In 2006, we initiated an extensive observational pro- 
gram to obtain multi-epoch high-angular resolution 
imaging observations of Galactic starburst clusters with 
the aim to study their internal dynamics and global mo- 
tions. Here, we present the results of our analysis of 
two epochs of HST observations of NYC separated by 
10.15yr. Accurate proper motions enable us to get a 
"clean" census of the cluster population by rejecting field 
stars and to study the internal cluster dynamics. 

2. OBSERVATIONS AND DATA ANALYSIS 

Two epochs of observations of NGC 3603 with the 
fiST/Wide-Field Planetary Camera 2 (WFPC2) sepa- 
rated by 10.15 yr are analyzed. We compare Planetary 
Camera (PCI) observations in F547M and F814W from 
epoch 1 (GO 6763, archival data) with our second epoch 
observations in F555W and F814W (GO 11193). With 
an image scale of 45. 5 mas pixel -1 , PCI provides the best 
point-spread function (PSF) sampling of the WFPC2 
cameras. While the first epoch observations were carried 
out in stare mode, for the second epoch we selected a 
four-point sub-pixel dither pattern to facilitate bad pixel 
recovery and to provide an improved PSF sampling. Ta- 
ble [1] summarizes the observations, including individual 
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Fig. 1. — Left: proper motions of stars in the NYC reference frame in Galactic coordinates. Right: histogram of the averaged one- 
dimensional proper motions of all stars (top), fitted by a two-component Gaussian (thick line), while the dashed line represents the 
distribution of the field stars. The lower panel depicts the observed proper motion dispersion as a function of stellar magnitude of the 
cluster member candidates. For 14.5 mag < mgi4 < 18.5 mag, we obtain a standard deviation o- b s lD = 184/^as/yr with an uncertainty of 
20/jas/yr (dashed line). 



TABLE 1 

Observing log 



Date 


Filter 




ttot(s) 


Stars 


30/07/97 


F547M 


3x1,12x10,8x30 


363 


772 


31/07/97 


F814W 


3x0.4,12x5,8x20 


221 


1163 


26/09/07 


F555W 


4x0.4,4x26,4x100 


506 


1048 


26/09/07 


F814W 


4x18,4x160 


712 


2014 



and total exposure times. 

Data reduction has been performed with IRAF/Pyraf. 
We combined the individual bias subtracted, flat- 
fielded images with identical exposu re times using 
multidrizzle (Koeke moer et al. 20021 ) . This corrects 
for velocity aberration and geometrical distortion includ- 
ing the 34th column anomaly based on the latest distor- 
tion correctior0__|. For the second epoch observations we 
applied 2x2 oversampling. 

Astrometry and photometry were derived from the 
drizzled images for each filter and exposure time setting 
using DAOPHOT (|Stetson 19901 ) with a Penny2 PSF 
varying linearly across the field. Near the faint end, 
where the photometric uncertainties start to increase, the 
star list was filled in by the results derived from the next 
longer exposure, and uncertainties assessed accordingly. 
The final number of detections in each band and epoch 
for a 5cr threshold above the background noise is listed 
in Table Q] Photometric calibration is based on the zero 
points provided by Dolphin (2000). Photometric correc- 
tion for charge transfer efficiency follows the recipe pro- 
vided by A. Dolphir0, and the astrometric correction is 
based on Equation 7 of Kozhurina-Platais et al. (2007) 
with the values for b\ , 62 , and &3 given in Table [2] (note 
the different pixel scale due to 2 x 2 oversampling for the 
second epoch). 

The combined imag e has distortions reduce d to a 
level of 0.02 pixel (see I Anderson fc King 2003bl ). Due 

3 http:/ /www. stsci.edu/hst/wfpc2/analysis/calfiles 

4 http:/ /ftp. stsci.edu/cdbs/uref/sadl946fu_idc. fits 

5 http:/ /purcell. as. arizona.edu/wfpc2_calib/ 



to the different orientation angle of 51° between the 
two epochs, we have to consider the uncertainty of the 
position induced by the residual geometric distortion, 
Cgco = .017 ± 0.001 pixel. For s inusoidal pixel phase er- 
rors (see lAnderson fc Ki ng 2000, their Figure 2), the sec- 
ond epoch dithering pattern with 0.5 pixel shifts largely 
cancels out the pixel phase error. As the first epoch 
was observed in stare mode, the pixel phase error has 
to be considered. With a typical amplitude of the sinu- 
soidal pixel phase error of ±0.02 pixel, the uncertainty 
to be included in our analysis amounts to an average 
residual uncertainty of cr pxp h = 0.01 3 ± 0.003 pi xel. Sim- 
ulations based on TinyTim PSFs (jKrist 19951 ) indicate 
that the positional PSF fitting uncertainty results in a 
centroiding error of erpsF = 0.013 ± 0.001 pixel. The 
effect of HST breathing on pixel scale was determined 
by measuring the separation of wide pairs of stars on 
frames obtained during different phases of HSTs orbit, 
resulting in CTbrcath = 0.009 ± 0.002 pixel. The com- 
bined contribution of these effects amounts to er crr = 

7geo + ^pxph + ^PSF + breath = 1-21 ± 0.18 mas. The 

observed proper motion dispersion has to be corrected 
for (Terr to derive the intrinsic velocity dispersion of the 
cluster members. 



TABLE 2 

Parameters of the astrometric correction 



Epoch 


£>i(px) 


6 2 (pxmag 1 ) 




1997.58 


0.094 


-0.0056 


6.17 x 10" 5 


2007.73 


0.122 


-0.011 


3.31 x 10~ 5 



As the orientation angles of the two epochs differ 
(A0 = 51°), the common field available to our analysis is 
a circle with a diameter of 30". A geometric transforma- 
tion based on a preliminary list of main-sequence (MS) 
cluster members is derived using IRAF/GEOMAP with a 
second-order polynomial. Stars detected in / band with 
photometric PSF fitting uncertainties cr p hot < 0.1 mag 
are matched after applying this transformation, and in- 
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dividual proper motions are calculated for each star. In 
this reference frame, cluster members center around (0,0) 
in the proper motion vector point diagram. The final 
proper motion table contains 829 matched stars. 

3. PROPER MOTIONS AND MEMBERSHIP 

In Figure [T] (left) we show the measured proper mo- 
tions for all stars with respect to the cluster reference 
frame in Galactic coordinates. The symmetrical distri- 
bution of all proper motions around (0,0) indicates the 
absence of any large relative motion of the cluster with 
respect to the field, i.e., the cluster follows the Galactic 
rotation curve. As NYC's Galactic longitude of 1=291.6° 
implies an almost tangential view into the Carina spi- 
ral arm, the distributions of cluster member and field 
star proper motions are superimposed. This proper mo- 
tion distribution is fitted by a two-component Gaussian 
with the wide function describing predominantly fore- 
ground stars, and hence non-cluster members (dashed 
line in Figure [T] top right), and the narrow component 
is interpreted as cluster member candidates. We calcu- 
late cluster membership probabilities P mem as described 
in Jones & Walker (1988) and consider stars as cluster 
members if P mem > 0.9. This significantly reduces the 
number of contaminating field stars, but due to the sim- 
ilar proper motions of cluster and field, some field stars 
might re main in the cluste r sample. Based on Besangon 
models (jRobin et al. 2 003). we estimate the number of 
field stars in our field of view to be 46 stars between 
16 mag < ni555 < 25 mag. Including stars with a mem- 
bership probability above 0.9 in our cluster sample, we 
subtracted a total of 58 stars as field stars. For the vari- 
able stars HST 12, 474, 481, and 574, studied by Moffat 
et al. (2004), we found membership probabilities P mem 
of 0.75, 0.98, 0.90, and 0.98, respectively, indicating that 
the latter three are likely cluster members. 

4. EXTINCTION, DISTANCE AND AGE 

The 777,555, 771555 — m-814 CMD is shown in Figure [2] 
(left). The distinction between candidate cluster mem- 
bers (small dots) and field stars (open circles) as de- 
scribed in Section[3]leads to a very well defined cluster se- 
quence with high-mass MS stars, intermediate-mass stars 
located in the pre-main-sequence (PMS)-MS transition 
region between 771555 = 18 mag and 20 mag, and more 
than 300 lower mass PMS stars down to 111555 — 24.5 mag 
(Figure [21 middle). The efficiency of the proper motion 
member selection is evident in the left panel of Figure 2, 
particularly among PMS and faint stars that could not 
be distinguished from cluster members from their colors 
alone. The foreground sequence blueward of the PMS 
does likely not belong to the cluster, suggesting a residual 
contamination of 18 stars with 20mag< 771555 < 24 mag. 

In the followin g, we assume solar metallicity 
(|Melena et al. 20*081) for the cluster and the relation be- 
tween absolute and selective extinction from Schlegel et 
al. (1998). The upper MS is we ll fitted by a 1 Myr Padova 
isochrone (Marigo et al. 20081 black solid line along the 
MS in the right panel of Figure [2|) for Ay =4. 7 mag and 
a distance modulus of 14.1 mag. For the analysis of the 
PMS-MS transiti on region and the lower mass PMS we 
use Siess models ()Siess et al. 20001) . transformed by Da 
Rio et al. (2009, S00/DR), as well as PISA-FRANEC 
models (Degl'Innocenti et al. 2008, PF09), transformed 



into the ob servational plane using A TLAS 9 model at- 
mospheres (|Castelli fc Kurucz 20 03). The best-fitting 
isochrones yield a distance between 6.6 kpc (PF09) and 
6.9 kpc (S00/DR) and a visual extinction Ay=4.7mag 
(PF09) and 4.6 mag (S00/DR), respectively for an age of 
1 Myr. We note that the derived selective extinction is in 
good agreement with E(B — V) = 1.25 mag as reported 
by Sung & Bessell (2004), though the absolute extinction 
value derived by us is slightly higher due to the use of 
the Schlegel et al. relations. The 1 Myr PF09 isochrone 
represents the PMS best, in particular at the PMS-MS 
transition region. 

At an age of 1 Myr, stars with masses between 3.5 and 
3.8 M© are expected to be i n the short-lived ra diative- 
convective (r-c) gap phase (|Mavne et al. 20071) . This 
phase corresponds to the formation of a radiative core 
in the interio r of the sta rs, due to the increasing central 
temperature |lbcn 1965J). We observe eight sources in 
the r-c gap at 18.5 mag< 777,555 < 19 mag (Figure and 
m 555 — ~ 2.25 mag (shown as diamonds in the right 
panel of Figure [2]). If their PMS nature is spectroscop- 
ically confirmed, this is the first identification of PMS 
stars in this interesting evolutionary stage. 

A previously unreported CMD feature is the appar- 
ent extension of the MS toward lower masses below 
the PMS-MS transition region (777,555 ^ 18.5 mag). 
Isochrone fitting to the MS turn-on yields an age of 
4 Myr. The derived age is consistent with recent esti- 
mates of the age of the two b lue supergiants Sher 23 and 
Sher25 (jMelena et al. 2 008). These stars might repre- 
sent an earlier epoch of star formatio n in the giant H II 
region (see also lSung fc Bessell 20041) . 

5. VELOCITY DISPERSION AND CLUSTER DYNAMICS 

The distribution <7 bsiD of the proper motions /i b s iD = 
Hi±m. ( as shown in the upper right panel of Figure [1] is 
a combination of internal velocity dispersion and instru- 
mental effects, resulting in cr bsiD = \J Vp^m + CT crr- I n 

the lower right panel of Figure Q] we show the observed 
one-dimensional proper motion dispersion as a function 
of stellar magnitude in bins of 1 mag for candidate clus- 
ter members. The velocity dispersion is constant for stars 
with 14.5 mag < mp 8 i4\v < 18.5 mag (~1.7-9M ). 

Correcting the observed one-dimensional proper mo- 
tion dispersion of cr bsiD = 184±20/ias/yr for the instru- 
mental effects discussed in Section[5]results in an intrinsic 
one-dimensional velocity dispersion £r pm iD = 141±27/ias 
yr _1 for stars brighter than 777,814 «18.5mag, assuming 
a negligible effect of binary orbital motions (Girard et 
al. 1989). This corresponds to ct c hd = 4.5 ± 0.8 km s _1 
at a distance of 6.75 kpc. The constant velocity disper- 
sion for stars in the mass range 1.7-9 M Q indicates a lack 
of equipartition of energy among cluster members. This 
provides a strong indication that NYC is far from virial 
equilibrium. 

Nevertheless, an upper limit of the cluster mass can 
be obtained by deriving the virial mass M nvn from the 
observed velocity dispersion (jSpitzer 19871 ): 



where 77 2.5 (weakly depending on cluster density 
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Fig. 2. — Left: 771555 —771814 vs. 771555 CMD of cluster (small dots) and field stars (open circles). Clearly visible are the MS and PMS loci 
at 771555 — 771814 ~1.5 and 3 mag, respectively. Center: cluster member candidates with P mom > 0.9. Apparent is the narrow PMS after 
the proper motion selection. Right: cluster member candidates with the best-fitting Padova i Myr (MS, black solid line) and PF09 1 Myr 
as well as S00/DR i Myr (PMS, dark and light gray solid lines, respectively) isochrones overplotted. Diamonds mark the stars in the r— c 
gap. The extension of the MS below the transition region is not covered by the younger isochrones, but is reproduced by a 4 Myr isochronc 
(PF09, dark gray dashed line; S00/DR, light gray dashed line). 



structure), rh is the half-mass radius, ct c 13d is the three- 
dimensional velocity dispersion, and G is the gravita- 
tional constant. 

NYC is mass segregated with its core ra- 
dius increasing with decrea sing stellar mass 
(jNiirnberger fc Petr-Gotzens 20021 ). A lower limit 
on rh for the high-mass stars as derived from HST data 
is comparable to the core radius of «0.2pc (Stolte 2003), 
whereas Harayama et al. (2008) based on the analysis 
of near-infrared adaptive optics data estimate rh = 
0.7-1.5 pc for stars in the mass range 0.5-2.5 M©. If 
we assume rh = 0.5 pc and a three-dimensional velocity 
dispersion of (j c \zt> = v3 x 4.5 ± 0.8km s _1 , we derive 
M dyn = 17600 ± 3800 M Q . 

Considering that this dynamical mass estimate pro- 
vides an upper limit, it is in agreement with photo- 
metric studies of NYC, which assigned masses to in- 
dividual stars, and estimated the total stellar mass 
to M c i g 10000 - 16000 M (|Stolte et al. 20061 : 
IHaravama et al. 20081) . 

6. SUMMARY AND CONCLUSIONS 

Based on two epochs of high-accuracy astrometric 
HST/WFPC2 observations separated by 10.15 yr, rela- 
tive proper motions of 829 stars were measured. A selec- 
tion of candidate cluster members with P mcm > 0.9 re- 
sults in a clean cluster CMD. The best-fitting isochrone 
yields an age of 1 Myr, a distance of 6.6-6.9 kpc, and 
Av=4.6-4.7mag for the PMS and intermediate-mass MS 
cluster members. 

Stars at the location of the short-lived radiative con- 
vective gap, which occurs at 3.5-3.8 M Q at the age of 
NYC, are identified for the first time. We find hints 
of a sparse young low-mass population with an age of 



~ 4 Myr, which might constitute an earlier generation 
of star formation in NGC 3603, and likely represents the 
low-mass counterparts to several blue supergiants in the 
vicinity of NYC. 

For the first time, the internal velocity dispersion of 
the starburst cluster NYC could be measured. For stars 
with masses 1.7 M© <M< 9M©, we determine a one- 
dimensional velocity dispersion of 141 ± 27/ias yr -1 , cor- 
responding to 4.5 ± 0.8 km s _1 at a distance of 6.75 kpc. 
From the fact that the velocity dispersion does not vary 
with stellar mass in this mass range, we deduce that 
NYC has not yet reached equipartition of energy. This 
is not entirely unexpected at the young age of the clus- 
ter, since its crossing time is estimated to be 1.4 Myr by 
Nurnberger & Petr-Gotzens (2002). 

The same might be true for many extragalactic star- 
burst clusters, where mass estimates rely on the mea- 
surements of velocity dispersions. If these clusters are 
also not yet in virial equilibrium, their masses might be 
systematically overestimated. Thus, NYC provides an 
important benchmark for our understanding of the early 
dynamical evolution and the long-term survival of young, 
massive stellar clusters in the Milky Way and in other 
galaxies. 
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